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ABSTRACT 

We study hydrogenated cations of two polycyclic hydrocarbon molecules as models of hy- 
drogenated organic species that form in the interstellar medium. Optical spectra of the hy- 
drogenated naphthalene cation H„-CioHg + for n = 1,2, and 10, as well as of the astro- 
biologically interesting hydrogenated proflavine cation H„-Ci3HnN3 + for n = 1 and 14, 
are calculated. The pseudopotential time dependent density functional theory is used. It is 
found that the fully hydrogenated proflavine cation Hi4-Ci3HnN3 + shows a broad spectrum 
where positions of individual lines are almost lost. The positions of lines, their shapes, and 
intensities change in hydronaphthalene and hydroproflavine cations showing that hydrogen 
additions induce substantially different optical spectra in a comparison to base polycyclic hy- 
drocarbon cations. One calculated line in the visible spectrum of Hio-CioHg + , and one in the 
visible spectrum of H-Ci3HnN3 + are close to the measured diffuse interstellar bands. We 
also present positions of near-UV lines. 
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1 INTRODUCTION 

Hydrogen, in its atomic and molecular form, is the most abun- 
dant of chemical species in the interstellar medium. The num- 
ber and size of detected interstellar molecules have been increas- 
ing last eighty years. Polycyclic aromatic hydrocarbons (PAHs) 
have been investigated as carriers of the diffuse interstellar bands 
(DIBs), because of unidentified infrared (UIR) emission features, 
and a s a source of an anomalous microwave emission (Tielens 
2008). It has been suggested that protonated PAH molecules 
play a catalytic role in t he process of the H2 formation in 
space (Bauschlicher 1998; lllirata et al.ll2004l) . Electronic transi- 
tions in the optical spectral region of protonated hydrogenated 
coronene, ovalene, pyrene, and circumpyrene have been recently 
calculated and discus s ed in the context of the DIBs p roblem 
dPafhak & Sarrel l2008t lHammonds. Pathak & Sarrel l2009t) . Spec- 
tra of isomers of protonated anthracene and phenanthrene have 
been measured in neon matr ices and studied by the time-dependent 
density functional theory (Garkusha et alj l201lh . We have se- 
lected here cations of naphthalene and proflavine to study the 
change of their overall optical spectra under additions of hydro- 
gen atoms. The prese nce of CipHs" 1 " in the i nterstellar medium 
has b een discussed (Iglesi as-Groth et a l. 2008; Galazutdin ov et all 
I20T 1). Several times in the past the evidence of the interstellar 
naph thalene cation was report ed and later shown to be prema- 
ture ((Galazutdinov et al. 2011). Proflavine is a substituted PAH 
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mol ecule with possible applications in astroph ysics and astrobiol- 
ogy dSarrell2006t feonaca & BilalbegovicfcOlOl) . 



A theoretical study of reactions has shown that the first H 
atom attaches to the naphthalene cation CioHg in exothermic 
reaction of about 60 kcal mol -1 , whereas the pr ocess for the 
second H atom is exothermic by 45 kcal mol -1 dBausch licher 
1998). Herbst and Le Page have found that the naphthalene 
cation efficiently associa tes with atomic hydrogen at all densities 
l lHerbst&Le Pagdll999h . Calculations have also shown that the 
addition of H at oms to the naphthalene cation pro ceeds with little 
to no barriers dRicca. Bakes & Bauschlichenl2007l) . Electronic ab- 
sorption spectra of H-CioHg + have been stu died by the photofr ag- 
ment spectroscopy and theoretical analysis dAlata et alTl20 1 OaH bh . 
It has been found that the protonated naphthalene cation absorbs 
in the visible part of the spectrum around 500 nm. A possibility 
that C10H8 and CiqHs + are carriers of DIBs has been studied 
( Salama & Allamandola 1 19921; Hirata, Lee & Head- Gordon! 1 19991; 
iKrelowski etal] 1200 ll : 



Malloci et all l2007bl) . The evidence of 



the naphthalene cation in the direction of the star Cernis 52 
in the Perseus molecula r cloud complex has been reported 
dlglesias-Groth et al.l2 008). The presence of some known DIBs has 
been confirmed in this study, and two new bands consistent with 
laboratory measurements on the naphthalene cation have been ob- 
served. It has also been proposed that hydrogen additions produce 
hydronaphthalene cations which contribute to the anomalo us mi- 
crowave emission i n this cloud ( Iglesias-G roth et al.l l2008). How- 
ever, a recent wo rk ( Galazu tdinov et al1l201 ll) has shown that this 
report (Igle sias-Groth et al.l2 008) is premature. The presence of the 
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naphthalene cation and related species in the interstellar medium 
deserves further studies. 

In addition to PAHs, the related molecules where C and H 
are subs tituted by other atoms are also studied in the astrophysica l 
context dHudgins. Bauschlicher & Alla mandola 2 005r , ISarrel 2006). 
Nitrogen is abundant in the interstellar medium. Its compounds 
have been detected in interstellar dust particles and meteorites. Or- 
ganic dye proflavine, C13H11N3 (3,6-diaminoacridine), has been 
proposed as one of molecules which act as a molecular "midwife" 
because of their property to accelerate DNA and RNA synthesis 
djain et al.l2 004). Therefore, the optical spectrum of proflavine and 
its ions is of astrobiological interest. Studies of such large biolog- 
ical molecules are impo rtant for understanding of prebiotic chem- 
istry dPuletti et alj201Cb . Organic dyes often exhibit strong lines in 
the visible spectrum. It has been shown that the optical spectrum of 
proflavine in aqueous solutions strongly depends on the pH value 
dDe Silvestri & Laportalll984l; Irlomem-de Mello et alj|2005) . Pro- 
tons can attach to nitrogen and carbon atoms in the molecule, and 
the state of proflavine in water and other solvents depends on pH 
of the solution. For example, it has been measured that the max- 
imum of absorption in the visible spectrum of proflavine (Amax) 
in water at room temperature, chan ges from 444 nm to 394 nm 
when pH changes from 7.0 to 14.0 dDe Silvestri & Laportalll984l) . 
The visible and UV optical spectra of proflavine and its ions have 
been recently studied using the p seudopotential time-dependen t 
density functional theory methods dBonaca & Bilalbeg ovicfcoiol) . 
Positions of spectral lines have been compared with DIBs, but with 
no definite conclusions. Because of the sensitivity of optical proper- 
ties of proflavine in water on the pH factor, we investigated here the 
impact of hydrogen additions on the optical spectrum of its cation. 

In this work we consider minimal and maximal hydrogena- 
tions of two organic cations. An addition of one, two and ten H 
atoms to the naphthalene cation, as well as one (in two positions) 
and fourteen hydrogen atoms to the proflavine cation are stud- 
ied. We use t he pseudopotential density functional theory (DFT) 
( Martin 2004) to determine the ground states of all base and hydro- 
genated species. Calculations of optical spectra of these systems 
within the pseudopoten tial time-dependent d ensity functional the- 
ory (TDDFT) methods dRunge & Gross! 1984 are described in sec- 
tion^ Results are presented and discussed in section|3] while final 
remarks are given in section [4] 



2 COMPUTATIONAL METHODS 



Optic al spectra are calculated using the Octopus code dCastro et alj 
2006). The geometries of all naphthalene and proflavine systems 
are minim ized independently b y the Quantum ESPRESSO DFT 
package (Giann ozzi et al J 12009). This approach is taken because 
of the fact that the Octopus is not in general suitable for a search 
of optimal geometries. Conditions as close as possible to simula- 
tions of the spectra within the Octopus code have been used in 
the Quantum ESPRESSO calculations (i.e., the local density ap- 
proximation (LDA) and corresponding pseudopotentials, see below 
and in Bonaca & Bilalbegovic 2010). The optimized geometries of 
cations are calculated without imposing any constraints, and they 
are used as inputs in the ground state and time-dependent calcula- 
tions in the Octopus code. 

Optical spectra are obtained using a time propagation method 
and the approximated enforced t ime-reversal symmetry algorithm 
dCastro. Marques & Ru bio 2004). The ground state is perturbed by 
an impulsive time-dependent potential and the time evolution is fol- 



lowed for 15.8 fs. In the time propagation TDDFT method used in 
this work, the width of sp ectral lines is inversely proportional to 
the total propagation time ()Castro et"ai . 2006; Marque s et al J2003l : 
iMalloci et al Jl2007bl : [Puletti et alJl201Cr) . The step of 0.0012 h/eV 
is applied. The absolute cross-section <j(lu) is obtained from dy- 
namical polarisability a{uf) which is calculated from the Fourier 
transform of the time-dependent dipole momentum of the system. 
Electronic spectra are calculated from 



<j(lj) — ^-Im a(w). 

TT 



(1) 



In this equation Im a(uj) is the imaginary part of the 
dyna mical polarisability. The Troullier-Martins pseudopoten- 
tials dTroullier & Martinsl LL993) and the TDLDA approxima- 
tion with the Perdew- Zunger exchange-correlation functional 
dPerdew & Zungerlll98lh are used, as well as the spacing of 0.13 
A in the real space method. The simulation cell is constructed 
from the spheres of the 4 A radii around atoms. The TDLDA ap- 
proximation in the pseudopotential TDDFT method shows a very 
good stability and produces results in an a greement with experi - 
ments, even for large biological molecules ( Marques et al.ll2003h . 
It has been found that the results of TDDFT calculations (within 
the Tamm-Dancoff approximation and using the BLYP functional) 
for the excitation ener gies of several PAH cations agree with exper- 
iments within 0.3 eV fairata. Lee & Head-Gordorl 1999b . We have 
found that calculations with the B3LYP functional (still in a devel- 
opment in the latest version of the Octopu s code) do not change 
subst antially the spectrum of proflavine dBonaca & B ilalbegovic 
l2010l) . It has also been found that spectra of proflavine obtained 
using the time propagation method in the Octopus code agree 
with results of the Lanczos chain TDDFT module in the Quant um 
ESPRESSO JWalker et alj|2006r . iBonaca & Bilalbegovic1l2010l) . It 
is known that the real-time propagation method of the Octo- 
pus code produces the whole optical spectrum up to the far-UV 
( Castr o et al. 2006). However, for energies above 10 eV only en- 
velopes of the spectra are accurate. In addition, there is little 
astronomical interest for a high-energy region of the spectrum. 
To facilitate a comparison with other studies of PAHs and or- 
ganic molecules o f astrophysical interest carried o u t by the same 
TDDFT method iMalloci. Joblin & Mulasl l2007al; IMalloci et alj 
l2007bl : lBonaca & Bilalbegovidl20ld : |Puletti et alj|2010h we show 
optical spectra up to 6 eV. 

The hydronaphthalene cation CioH^ and the dihydronaphtha- 
lene cation CioH^,, with additional H atoms attached at the same 
positions as in Bauschlicher (1998), are studied. Two situations of 
a minimal hydrogenation, where one hydrogen atom is added to the 
proflavine cation, are also investigated here: an additional H atom 
attached as the second hydrogen atom at the central carbon atom 
(labelled as (I)), and at the opposite nitrogen atom in the same 
central ring (labelled as (II)). We also study fully hydrogenated 
species of naphthalene and proflavine cations: Hio-CioHs + , and 
Hi4-Ci Hi3N 3 + . 



3 RESULTS AND DISCUSSION 

The optimized structures of all investigated species of naphtha- 
lene and proflavine are shown as insets in Figs. 1-3. In contrast 
to the naphthalene and proflavine molecules, carbon skeletons of 
fully hydrogenated Hio-CioHs + and Hi4-Ci3HnN3 + (shown in 
Fig.|2jc) and Fig.[3jc)) are not planar. In addition, hydrogen atoms 
attached in pairs on peripheral carbon atoms are positioned below 
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Figure 1. Optical spectra of the naphthalene cation: (a) after the geometry 
optimization without constraints, (b) for the cation with the optimized struc- 
ture of the neutral naphthalene. Insets: Black balls represent carbon atoms, 
whereas white ones mod el hydrogen a toms. Structures are visualized using 
the XCrySDen program (Kokali 2003). 
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Figure 2. Optical spectra of hydrogenated species of naphthalene cation: 
(a) H-Ci H 8 + , (b) H 2 -Ci H 8 + , (c) H 10 -Ci H 8 + . 



and above the middle plane of carbon rings. Average C-C distances 
increase from 1.39 A in the optimized naphthalene cation to 1.51 
A in Hio-CioH 8 + . Average N-C distances increase from 1.34 A in 
Ci 3 HnN 3 + to 1.41 A in Hi 4 -Ci 3 HiiN 3 + . Average C-C distances 
increase from 1.40 A in the proflavine cation to 1.51 A in H14- 
Ci 3 HnN 3 + . Hydrogen atoms are also attached on the C atoms 
in pairs above and below rings in optimized H-CioH8 + (one pair, 
inset in Fig. |2ja)) and H2-CioHg + (two pairs, inset in Fig.|2jb)). 
However, carbon skeletons are planar in both these structures with 
a minimal hydrogenation, as well as in H-Ci 3 HnN 3 + (in both po- 
sitions I and II, insets in Fig.[3ja) and Fig.[3lb)). The similar con- 
clusion about a planarity of H-Ci oHs + has been reached by Alata 
and coworkers dAlata et al.l2010al) . 

Optical spectra presented in this work are calculated start- 
ing from structures optimized separately for all pure and hy- 
drogenated cations, and without any constraints in geometri- 
cal optimizations. We found that the fully optimized naph- 
thalene cation is still a planar structure. However, the charge 
and small displacements from atom positions of the neutral 
naphthalene produce small differences in the spectrum. Opti- 
cal spectra for the fully optimized cation and the structure 
where one electron is remo ved from the optimized geome- 
try o f a neutral naphthalene IINiederalt. Grimme & P everimhoff 
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1 19951 : Iffirata. Lee & Head-Gordonl 1 19991 ; iMalloci et al j 12007b) 
are shown in Fig. [TJ The geometry used as an input for 
Fig. [ji b) is taken from the Theo retical spectral database of PAHs 
(Malloci, Joblin & Mula3 l2007al) . Our result for the optical spec- 
trum of the naphthalene c ation in this situation agrees w ith one 
presented in the database (Mall oci. Joblin & Mu las 2007a). Lines 
for spectra shown in Fig.[TJare compared in Table [TJ 
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Figure 3. Optical spectra of hydrogenated species of proflavine cation: 
(a) H-Ci 3 HnN 3 + , additional H atom is attached at the central C atom 
(position I), (b) H-Ci 3 HnN 3 + , additional H atom attached at the central 
N atom (position II), (c) Hi4-Ci 3 HnN 3 + . Insets: The nitrogen atoms are 
labelled by "N" on the right-hand side of corresponding balls. 
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Optical spectra of hydrogenated protonated naphthalene 
species are shown in Fig. [2] Peaks of the naphthalene cation move 
in H-CioH8 + , H2-CioH8 + , and Hio-CioH8 + , and their intensi- 
ties change. Optical spectra of hydrogenated protonated proflavine 
species are shown in Fig. [3] The most notable fact is that the spec- 
trum of H14-C13H11N3 4 " (Fig.[3jc)) is broad and this broadening is 
much more pronounced than in the fully hydrogenated naphthalene 
(Fig.flc)). 

Spectral lines of all naphthalene and proflavine related sys- 
tems are presented in in Table[TJand Table|2] Visible spectra of hy- 
drogenated naphthalene cations extend up to ~ 500 nm, whereas 
long tails of visible spectra of hydrogenated proflavine cations 
spread above 800 nm. The photofragmentation spectroscopy mea- 
surements by Alata and coworkers have also been shown that 
H-Ciq H8 + absorbs in the visible, around 500 nm dAlata et all 
l2010al lbh. We also present near-UV lines (above 200 nm) to facil- 
itate a comparison with corr esponding spectral mea surements of 
interstellar organic materials dKwok & Sanfordl2008l) . Such exper- 
iments are designed to be carried out on the Cosmi c Origins Spec- 
trograph (COS) of the Hubble Space Telescope dOsterman et alj 
|2011|) . The near-UV channel of COS shows a high sensitivity in the 
region between 200 nm and 300 nm where lines of hydrogenated 
protonated naphthalene and proflavine exist. The lines of 28 1 nm 
and 422 nm in the spectrum of the proflavine cation [Bonaca2010] 
move to the 273 nm and 454 nm and three new lines appeared 
between these two in the spectrum of H-Ci3HnN3 + (position I). 
Rather strong lines exist at 201 nm, 222 nm, and 244 nm. An addi- 
tion of the H atom to the central N atom produces lines between 206 
nm and 470 nm, as shown in Table [2] The change of optical lines 
is the most obvious in the spectrum of Hi4-Ci3HnN3 + , where the 
UV lines exist at 216 nm and 253 nm. Then intensity drops, and the 
lines at 352 nm, 420 nm, and 496 nm are very weak. 

It is known that the majority of DIBs are present in the 
atomic hydrogen gas dHerbigl [l993l: ISnow & McCalll |2006|) . In 
Table [3] we com pare positions of several calculated lines wit h 
the closest DIBs dHobbs et alj|2008l , 120091 : iTuairisg et alj|2000h . 
The best agreement is for a line of the fully hydrogenated 
naphthalene cation HioCioH8 + , and a line of H-C13H11N3 4 " 
(position II). Alth ough TDDFT calculations show only trends 
in DIBs position s (Malloci . Joblin & M ulas 2007a; Mallo ci et"al] 
l2007bl : IPathak & Sarrell2008l ; iHammonds. Pathak & Sarrell2009l) . 
the agreement presented in Table [5] is good and deserves a fur- 
ther experimental investigation. In the recent report of the high 
resolution spectroscopy observation of CioHs 4 " in the Cernis 52, 
the importance of hydrona phthalene cations has been pointed out 
(Iglesias-G roth et afll2008h. Although this r eport has been shown 
to be premature dGalazutdinov et all 1201 lh . the presence of the 
naphthalene cation and its hydrogenated forms in the interstellar 
medium remains an important subject for further investigations. 



4 CONCLUSIONS 

Studies of spectra of polycyclic aromatic hydrocarbon cations with 
attached hydrogen atoms are important for the DIB and UIR prob- 
lems, an anomalous microwave emission, as well as for general 
properties of organic matter in the interstellar medium. We study 
the changes in electronic absorption spectra induced by hydrogen 
additions to the naphthalene and proflavine cations using the time- 
dependent density functional method in its pseudopotential version. 
Calculated spectra are based on the overall density of electronic 
transitions and show that hydrogen additions substantially change 



intensities, shapes, and positions of optical lines in a compari- 
son with spectra of base cations. Therefore, additions of hydrogen 
atoms is important in astrophysical applications of optical spectra 
of organic molecules. Similar conclusions have been found for pro- 
tonated hydrogenated coronene, ovalene, pyrene, and circu mpyrene 
dPathak & Sarre 20081; IHammonds. Pathak & Sarrel 12009). Calcu- 
lated lines of protonated hydrogenated naphthalene and proflavine 
are compared with measured DIBs. An especially good agreement 
exists for the visible line of the fully hydrogenated naphthalene 
cation and two visible lines of the hydroproflavine cations. The 
lines in the near-UV spectral region are also presented. Our calcula- 
tions should give guidelines for the change of near-UV and visible 
spectral lines for similar larger cations and their derivatives under 
the process of hydrogenation. 
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Table 1. Optical spectral lines of the naphthalene cation and hydrogenated naphthalene cation species. The star labells data for the naphthalene cation with the 
optimized geometry of the neutral naphthalene molecule, whereas all other data are obtained for optimized cations. Because of a limited accuracy of TDDFT 
method, the calculated results are presented rounded to whole numbers. 
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Table 2. Optical spectral lines of the proflavine cation and hydrogenated proflavine cation species. Labels (I) and (II) are for H atoms attached, respectively to 
the central carbon and central nitrogen atom. Because of a limited accuracy of TDDFT method, the calculated results are presented rounded to whole numbers. 
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Table 3. Lines calculated using the time-dependent density functional theory (TDDFT) method compared to the nearest DIBs. Because of a limited accuracy 
of TDDFT method, the calculated results are presented rounded to whole numbers. FWHM of experimental lines are also shown. 
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